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Abstract 

We consider a three-dimensional Ising model in a transverse magnetic field, 
h and a bulk field H. An interface is introduced by an appropriate choice of 
boundary conditions. At the point {H = 0, h = 0) spin configurations corre- 
sponding to different positions of the interface are degenerate. By studying 
the phase diagram near this multiphase point using quantum-mechanical per- 
turbation theory we show that that quantum fluctuations, controlled by h, 
split the multiphase degeneracy giving rise to an infinite sequence of layering 
transitions. 

PACS numbers: 75.10.-b, 75.50.Rr, 68.45. Gd 
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I. INTRODUCTION 



There is a considerable body of literature discussing the way in which interfaces depin 
from surfaces[ |l|. Of particular interest to us here is the situation below the roughening 
transition when the interface is smooth and can depin from the surface through a series of 
first-order layering transitions. 

This possibility was first pointed out by De Oliveira and Griffiths [ [J for a model in 
which the layering is driven by the competition between a long-range bulk interaction and 
a surface field. Here the layering transitions exist even at zero temperature. Later Duxbury 
and Yeomansf [|| showed that, if the position of the interface relative to the surface was 
degenerate at zero temperature, the degeneracy could be split by thermal fluctuations giving 
an infinite series of layering transitions at finite temperatures. The stable interface position is 
determined as a balance between the binding effect of a bulk field and the entropic advantage 
for the interface lying further from the surface. 

The transverse Ising model [ |J was first introduced by de Gennes [ || in connection 
with ferro-electric materials. Recently Henkel et. al.[ || have discussed the behaviour of 
a domain wall in this system. The interesting questions concern the effect of quantum 
fluctuations, mediated by the transverse field, on the behaviour of the interface. Their work 
considers one dimension where the interface is rough. Very different behaviour is likely for a 
smooth interface. Therefore in this paper we consider the behaviour of the three-dimensional 
transverse Ising model below the roughening temperature. We find that a zero-temperature 
multiphase point can be split by quantum fluctuations and that, for a non-zero transverse 
field, there are an infinite number of stable positions for the interface as a bulk field passes 
through zero. 

The next section of the paper defines the model and gives a qualitative discussion of 
its properties and our approach. Quantitative details of the calculation, which is based 
on quantum mechanical perturbation theory, are given in Sec. III. Our conclusions are 
summarized in Sec. IV. 
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II. QUALITATIVE REMARKS 



The Hamiltonian of the three-dimensional transverse Ising model we shall consider is 

n=-j J2 T, s *(hj)Sz(hf)-j E J2s z (i,j)s z (i + i,j) 

i=l,L(j,j>) i=0,L j 

- E Y,\hS x (i,j) + HS z (i,j)} -Kj2(Sz(0,j)-S z (L + l,j)) , (2.1) 

i=l,L j 1 J j V 1 

where i labels two-dimensional planes and nearest neighbours within a plane. The 

parameter K — > oo is used to impose appropriate boundary conditions, namely to fix the 
spins at the surface (i = 0) to be up and those in the last layer (i — L + 1) to be down. 
These boundary conditions will create a domain wall, or interface, in the system separating 
layers of up and down spins (see Fig. [j]). Our aim is to construct the phase diagram which 
gives the position, k, of the interface (defined as in Fig. ^|a.) as a function of the uniform 
field H and the transverse field h. We shall consider the limit L — > oo and zero temperature 
where the interface is flat. 

As a first step in understanding the phase diagram we consider the situation for h = 0. 
In this case it is clear that for positive H, k = oo whereas for negative H, k = 0. We shall 
call these phases, which are illustrated in Fig. |^, R and L respectively. For if =0 the energy 
is independent of k, so that all interface positions are degenerate. It is known that such a 
degeneracy can be lifted by either thermal fluctuations [ 0] or by quantum fluctuations[ ||. 
In more general contexts this removal of degeneracy has been referred to as ground state 
selection [ |9| following the work of Villain and Gordon [ |TT| and Shenderf [10| . Here we will 
consider the effect of quantum fluctuations. 

There are two possibilities: in the first, quantum fluctuations due to nonzero h cause 
the transition from L to R to be discontinuous with no intermediate states; in the second 
this transition occurs through a sequence of intermediate states in which k increases mono- 
tonically. This sequence can be finite, so that there is a first-order transition from a state 
with k = k ma-x to R, or it can be infinite. In general, one expects the first possibility when 
the effective interaction between the surface and the interface is attractive and the second 



when this effective interaction is repulsive. Our results indicate that quantum effects give 
rise to the second possibility and that the sequence of layering transitions is probably an 
infinite one. It is easy to see that even when h is nonzero, the R phase is stable whenever 
H is positive. As we will show, the stability of the L phase requires that H < —Ch 2 j/Jl, 
where C is a constant. 

To determine the interface position we must calculate the ground state energy as a 
function of H . To do this we assume perfectly flat interfaces and apply perturbation theory 
to calculate the energy of the state when the interface is at position k. If we were dealing with 
a finite system, then perturbation theory would introduce coupling in finite order between 
states when the interface is at different positions. However, this tunneling effect disappears 
in the thermodynamic limit. Our calculations are valid for h ^ J. We also impose the 
condition j <C J , although as long as we remain in the regime where the interfaces are flat, 
this restriction is probably inessential. 

The energy per layer spin e(k) of the state when the interface is at position k can be 
written 

e(k) = e (H) - kH + E k (H) , (2.2) 

where eo(H) is the energy per layer spin for the k = phase. The most important dependence 
on H is in the term —kH . H also contributes to the energy denominators that appear in 
perturbation theory, which are typically of the form 2 J + j + H. However perturbative 
contributions to E k (H), when expanded in powers of H, lead to corrections which are of 
relative order Hj J or smaller. But since we will only be interested in H in the range 
—Ch 2 /J < H < 0, these corrections are smaller than of relative order {h 2 / Jq) <C 1, which 
we may ignore. Therefore, in Eq. ( |2.2| ), we may evaluate E k (H) at H = 0. We shall find 
that E k is a positive and decreasing function of k. This result leads to an infinite sequence 
of phase transitions. The critical field separating the phase k from k + 1 is given by 

— HI = E k — E k+1 = AE k , (2.3) 

where, as argued above, the leading order calculation of E k can be obtained for H = 0. 
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III. CALCULATION OF AE K 



From Eq. ( f2.3| ) it is clear that we only need to keep track of terms in the ground-state 
energy which depend on k. In other words we need to ascertain how the corrections to the 
ground-state energy which are perturbative in h depend on the location of the interface. 
For convenience we now transform to occupation number operators. For spins that are up 
(down) we write S z (i,j) = 1/2 - n i:j (S z (i,j) = -1/2 + n itj ). Also S x (i,j) = (a\ d + a itj )/2, 
where the operator a\j (a it j) creates (destroys) a Bose excitation at site % in the jth layer 
and riij = a it ja\j. Strictly speaking we should not allow more than one excitation to exist 
on a single site. To enforce this restriction we include a term of the form A J2i. 
where A — > oo. Normally, it is difficult to take full account of such a hard-core interaction. As 
we will see, we accommodate this constraint by never involving matrix elements connecting 
to a state in which there is more than one excitation at any site. Therefore, setting H = 0, 
we are lead to the following bosonic Hamiltonian, when the interface is at position k: 

OO , s oo oo 

H B = E + ( 2J o n i,j - (V 2 ) [a{j + a i,j] -EE Jorkjriiji + A^E a lj a {j a ij a i,j 

i=i j v J i=\ (j,j>) i=\ j 

oo 

-j EW 1 / 2 )^ + + n iJ n i+i,j} + J ^2l~ n k,j ~ n k+ij + 2n kJ n k+1>j ] , (3.1) 

i=0 j j 

where E is the unperturbed energy of the k = phase, A — > oo, and for K — > oo, we may 
set noj = 0. We write this Hamiltonian as 

H B = E + Ho + V 1 + V 2 + V 3 + V, , (3.2) 

where 

^o = A EE n iJ ' (3-3) 

* 3 

with A = 2 J + j, 

Vi = - EEiv+ in «.i n »+i,i _ E EZ n iJ n hj' (3- 4 ) 

^ = -(V2)EEli + ^] . (3-5) 

* 3 



V 3 = -J Y.( U kJ + n k+l,j) , 



(3.6) 



j 



where 



j for i 7^ k 



(3.7) 



Ji,i+1 



j for i — k 



and (with A — > oo) 



^ = aEE 



@ J i,j@ J i,j@ J i,j®'i,j • 



(3.8) 



* 3 



We now consider how the perturbative contributions to the energy depend on the various 
coupling constants. To carry out this discussion it is convenient to introduce a diagram- 
matic representation of the contributions to the perturbation expansion. Each term of V\ 
proportional to j rs n r n s , where r (and similarly s) denotes a position label of the form i,j, 
is represented by a line joining the two interacting sites r and s and j rs = Jq or j rs = j 
depending on whether sites r and s are in the same plane or are in adjacent planes. The 
perturbation in V 2 proportional to a r (aj) is represented by a minus (plus) sign at the site 
r. However, for simplicity, since each site involved with any of the preceding interactions 
must be excited (i. e. must have both a "+" and a "-" associated with it), we have not 
explicitly shown "+"'s and "-"'s in Fig. |||. The term in the perturbation V3 proportional to 
n k,j ( n k+i,j) is represented by a circle attached to the site k,j (k + 1, j). Any term in per- 
turbation theory which does not involve V4 can be constructed from these elements. Some 
simple examples are shown in Fig. [|. 

We now define what we mean by connected terms. Any term which involves only a single 
site is connected. Terms which involve more than one site are connected only if all such 
sites are connected with respect to lines representing terms of V\. If this is not the case, the 
term will be called disconnected. Thus diagrams (a) and (d) of Fig. |3] are disconnected. 

We now establish that contributions from disconnected configurations of lines vanish. 

Consider a disconnected diagram V which consists of two disjoint components, r a and 
r#. The contribution of of this diagram is unchanged if we were to treat perturbatively the 
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system S(Ta + Tb) in which all coupling constants j rs not in or Tb are set to zero. But 
because Ta and T B are disjoint systems, we have E(Ta + T b ) = E(Ta) + E{Tb)- This result 
indicates that there are NO disconnected terms in the ground state energy which involve 
simultaneously an exchange constant from one component Ta and an exchange constant 
from the other component Tb- Thus disconnected diagrams can be omitted from further 
consideration. 



To evaluate AEk it is apparent from the form of Eq. ( |2.2|) that we only need to keep 
contributions which appear when the interface is at position k + 1 but NOT when it is at 
position k. Note that if the diagram does not involve the interface potential V3, it has no 
dependence on position and will make equal contributions in the two cases. So AE^. has 
contributions from diagrams which both a) involve V3, the interface potential, and b) can 
occur when the interface is at position k + 1 but not at position k. Such a connected diagram 
must involve sites in rows 1, 2, ... k + 1. This can be done with diagrams involving k or 
more lines. As we will see in a moment, the dominant contribution involves the least number 
of lines. To take the least number of lines, means that we take k lines representing j and 
no lines representing J . Thus the dominant diagrams are linear chain diagrams, and we 
may therefore consider a Hamiltonian in which the index j in Eqs. ( |3.3| )- P^ ) is omitted. 
We illustrate the diagrams with the minimum number of lines which contribute to AEk for 
k = and k = 1 in Figs. ^ and [5] respectively. 

To check that it is indeed the linear chain diagrams which give the lowest order contribu- 
tion to AEj, it is necessary to give a more detailed analysis of the contribution of a diagram 
involving, say, p different lines representing V\ and, as we explained, necessarily involving at 
least one interface potential term V3. Note that both these perturbations, V\ and V3, involve 
occupation numbers n r , which vanish when there is no excitation at site r. For every site 
r involved in a V\ or V3 interaction, it is necessary to create an excitation so that n r can 
be evaluated in a virtual state in which n r — 1. Subsequently, in order to get back into the 
ground state we must destroy the excitation on the site r. Thus in all a diagram involving 
p different lines will involve p+1 sites and therefore give rise to a perturbative contribution 



to the energy which is of order S P E, where 

5 P E = h ( - 2p+2 ^f 1 J P2 (j/A)/A 2p+1+Pl+P2 , (3.9) 

where of the p lines, p\ are associated with j and P2 with Jo [see Eq. ( |3.4j )1. In writing 
this equation we included the factor (J / A) to take account of the necessary factor of V3. 
At this point, it is clear that to have a diagram which occurs for position k + 1 but not 
for k it is best to invoke a linear diagram, and not one which reaches more than one row 
perpendicular to the interface. Unnecessary factors of Jo in V\ will give rise to additional 
factor of Joh 2 /A 3 <C 1. So, we conclude that to leading order 

AE k = C k J 2 [jhl/J*} k+1 , (3.10) 

where C k is a constant which must be determined by an explicit calculation and to leading 
order we set A = 2J . 

Now let us carry out a detailed calculation for the simplest case, namely for k — 0. From 
Fig. H] we obtain 

AE = -E(A1) , (3.11) 

where E(A1) is the contribution to the energy in diagram A\ of Fig. |4|. Thus[ |12| 

AE = -(0|F 2 i^y 2 |0) 

= -(0\(-h/2)a 1 ^(-ja\a 1 )±(-h/2)a\\0) 

= (h 2 /A)(j/A 2 )=jh 2 /(16J 2 ) . (3.12) 

Here and below the excitation energies £ will be — rA = — 2r J , where r is the number of 
excitations in the virtual state. 

Next we calculate AJ^i = —E(E2) from diagram E2 of Fig. 0. Here we have to sum 
over the different orderings of the perturbations. 

AE 1 = 

8 



-<o| 



1 



(-h/2) ai ~{-h/2)a2 + {-h/2)a 2 j{-h/2) ai 
-h/2)a\j(-h/2)a\ + {-h/2)a\^(-h/2)a\ 



-jnin 2 )-(-jn 2 ) + {-jn^-^-jn^) 



10} 



■<o| 



{-h/2)a^(-h/2)a 2 + {-h/2)a 2 j{-h/2) ai 



^Hnin 2 )^(-h/2)a\^(-jn 2 )^(-h/2)al\0) 



-(0|(-V2)a 2 |(-jn 2 )|(-V2)ai|(-jnin 2 )i 



(-h/2)a\^(-h/2)4 + (-h/2)4^(-h/2)a\ 



|0>. 
(3.13) 



When simplified this yields 

= 2(V2) 4 J 2 /A 5 = (1/256)^/ J 5 . (3.14) 

This calculation is hard to extend to AE^ for larger k using naive perturbation theory. 
We need a more powerful formalism, namely Matsubara diagrams[ In this formalism 
one has diagrams constructed from the following elements. The perturbations V\, V 2 , and 
V3 are represented by vertices as shown in Fig. ^. Each such vertex carries the appropriate 
factor ((— h/2), —jS r) k, and —j, respectively, where 8 r ^ is the Kronecker delta. Also note 
that at leading order j^i+x = j, because we never invoke the term jk,k+i-)- Lines labeled with 
the same index are joined and a sum is taken over all topologically inequivalent connected 
diagrams. Each line represents a Green's function (z u — A) -1 . All the indices are summed 
over. The #'s are summed over the Matsubara frequencies z v = 2viti/ (ksT) where v runs 
over all integers positive and negative. One enforces conservation of z, that is for each 
vertex the sum of all incoming z's minus the sum of all outgoing z's must equal zero. For 
the present case, this conservation law mean that at any vertex which has only one line 
entering or leaving the corresponding z must be zero. One can quickly see that the z's for 
all lines have to be zero. So, in fact, there is no sum over z to be done. (Such sums normally 
lead to the Bose occupation number factors.) 

From Fig. [7] we see that 

A£ = -E a , (3.15) 
where E a is the energy from diagram (a) of Fig. [7]. (A similar notation for the other diagrams 



is used below.) Thus 

AE = _(-V2)(-j)(-V2)(-A)- 2 = h 2 j/{AA 2 ) = h 2 j/(16J 2 ) , (3.16) 

as before. To obtain this result, note that diagram (a) has two filled circle vertices (each 
carrying a factor —h/2), one triangle (carrying a factor — j), and two lines, each of which 
carries a factor (z — -Ej) 1 = A -1 . Also from Fig. |7| we have 

AE 1 = -{E b + E e ) = -2E b . (3.17) 
Since diagram (b) has four h vertices, two j vertices, and five lines, we have that 

AEl = _ 2 (_/ i /2) 4 (-i)(-i)(-A)- 5 = /i 4 j 2 /(8A 5 ) = h 4 j 2 /(256J 5 ) . (3.18) 
Finally, from Fig. [7| we have 

AE 2 = -{E d + E e + E f + E g ) = -AE d . (3.19) 
Since diagram (d) has six h vertices, three j vertices, and eight lines 

AE2 = _4(_/ i / 2 )6(_ i )(_ i )2 ( _ A) -8 = h *jS/(i 6A *) = ^-3/(4096^) • (3-20) 

Evidently, the general result is 

AE k = U 3 h 2 /lQJl) k+1 . (3.21) 

IV. CONCLUSIONS 

Equation ( |3.21| ) indicates that the boundary between phases with k = p and k = p + 1 
is given to leading order by H* = — Jo(j h 2 /16Jq) p+1 . The resulting phase diagram is shown 
schematically in Fig. |8|. 

The analysis presented above was based on retaining only the leading-order (in H/Jq and 
j / Jo) term in the surface-interface interaction. Although we cannot rule out the possibility 
that the neglected higher-order interactions could become dominant for very large k (for fixed 
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h and Jo), we do not expect to observe any qualitative corrections to the phase diagram in 
this limit. This is because there are no competing interactions which would make correlation 
functions oscillatory at large distance and therefore it seems implausible that the positive 
sign of AE p can be changed by the neglected higher-order terms[ |14|| . 

To summarize we have shown that quantum fluctuations can stabilize an infinite sequence 
of layering transitions in a three-dimensional transverse Ising model. 
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FIG. 2. a) Configuration of the spins when the interface is at position k. Layer is constrained 
to be up and the layer at infinity is constrained to be down. In panels (b) and (c) we have 
represented the configurations R and L. 
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FIG. 3. Representation of perturbative contributions to the ground state energy. The dashed 
line corresponds to the location of the interface, (a) A contribution which cannot occur, (b) A 
term in which a given j is taken to second order, (c) The contribution to lowest order in j which 
involves four sites. (d),(e),(f) Terms which involved factors of h 2 . (f) A term which involves the 
Jo interaction. 
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FIG. 4. Perturbative contributions to the ground-state energy used to evaluate AEq = E$ — E\. 
Panel AO corresponds to the case when the interface is at k = 0; panels Al and Bl to the case 
k = 1. 
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k = 2 

FIG. 5. Perturbative contributions to the ground-state energy used to evaluate AE\ = E\ — E2. 
Panels Al - Dl are for the case when the interface is at k = 1. Panels A2 - E2 correspond to 
k = 2. The dashed line represents the position of the interface. Contributions from diagrams Al, 
Bl, CI, and Dl are equal to those from diagrams A2, B2, C2, and D2, respectively. So AEi is 
determined by E2. 
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FIG. 6. Vertices for Matsubara diagrams: (a) transverse field vertex —{h/2){a\ n ) (top) and 
-(h/2)a m (bottom); (b) -j5 r)k a\a r ; (c) -jrijUj+i = ja< j a i j+1 a j+ ia j . 
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FIG. 7. Matsubara diagrams. The vertices are explained in the preceding figure. The lines are 
labeled by an index i to represent the Green's function (z — A) -1 . This calculation is for a linear 
system with sites labeled 1, 2, .... These diagram give — AE m following the reasoning of Figs. 4 
and 5. Diagram (a) is for m = 0, diagrams (b) and (c) are for m = 1 and diagrams (d)-(g) are for 
m = 2. No new topology is obtained by reversing the direction of the line labeled "1". However, 
one may independently reverse the direction of all other lines giving rise to a degeneracy 2 m . 




FIG. 8. Schematic representation of the phase diagram for interface unbinding transitions in 
the transverse Ising model. 
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